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CHAPTER  I 


INTRODUCTION 

Many  of  the  materials  in  common  usage  today  are  orthotropic 
in  their  elastic  properties.  The  different  types  of  wood  are  the  most 
familiar;  another  less  familiar  is  rolled  aluminum  in  which  the  rolling 
process  seems  to  alter  the  elastic  properties  of  the  aluminum  and  produce 
one  modulus  of  elasticity  in  the  direction  of  the  rolling  and  another 
perpendicular  to  the  direction  of  rolling.  As  an  illustration  of  the 
properties  of  an  orthotropic  material  consider  Sitka  spruce:  the  modulus 
of  elasticity  parallel  to  the  grain  is  1,679  x lC^  lbs.  per  square  inch 
while  the  modulus  of  elasticity  perpendicular  to  the  grain  is  0.076  x 
10^  lbs.  per  square  inch.  It  is  evident  that  serious  discrepancies 
could  arise  if  the  orthogonality  of  the  elastic  properties  of  a material 
such  as  Sitka  spruce  were  not  taken  into  account  when  it  is  put  to  some 
practical  use.  This  is  especially  true  when  a material  is  to  be  used 
where  it  will  be  stressed  close  to  its  elastic  limit  such  as  in  an  aircraft. 
Due  to  the  need  to  eliminate  all  excess  weight,  a member  is  designed  so 
that  the  material  is  used  to  its  utmost  advantage.  It  might  be  added  that 
quite  often  for  some  materials  the  limiting  stress  is  not  determined  by 
the  modulus  of  elasticity  but  by  the  modulus  of  rigidity,  that  is,  the 
shear  stress  reaches  its  limiting  value  before  the  normal  stresses  reach 
their  limiting  values. 

Consider  then  a semi-infinite  plate  oriented  so  that  one  of 

the  axes  of  elastic  symmetry  is  parallel  to  the  boundary  of  the  plate. 

Such  a plate  might  be  a large  piece  of  plywood.  Let  the  axis  along  the 
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boundary  of  the  plate  be  the  x-axis  and  let  the  y-axis  extend  in  the 
positive  direction  into  the  plate.  The  plate  will  be  subjected  to  a 
load  which  will  act  in  the  plane  of  the  plate  and  perpendicular  to  the 
boundary  of  the  plate.  The  resulting  condition  will  be  treated  as 
plane  stress. 

The  components  of  stress  and  strain  in  an  orthotropic  plate 
in  a state  of  plane  stress  are  connected  by  the  following  relations:1 
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where  the  x-  and  y-axes  are  the  axes  of  elastic  symmetry.  In  these 
equations  Ex  and  Ey  are  Young's  moduli  of  elasticity  in  the  x-  and  y- 
directions  respectively.  Poisson's  ratio  (T'xy  is  the  ratio  of  the 
contraction  parallel  to  the  y-axis  to  the  extension  parallel  to  the 
x-axis  associated  with  a tension  parallel  to  the  x-axis  and  similarly 
for  C yx.  The  quantity  y^fxy  is  the  modulus  of  rigidity  associated 
with  the  directions  of  x and  y.  It  should  be  mentioned  that  these 

1H.  W.  March,  Stress-Strain  Relations  in  Wood  and  Plywood 
Considered  as  Orthotropic  Materials,  Porest  Products  Laboratory  Report 
M>VRl503,  february,  TOi,  p.  2. 
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relations  hold  only  for  small  deformations  where  all  of  the  stresses 
produced  are  within  the  elastic  limits  of  the  material. 

The  equations  of  equilibrium  for  plane  stress  are 


•J&-  x . ■+  — V -o 


These  equations  are  deduced  from  the  conditions  of  static  equilibrium 
with  no  reference  to  the  relations  between  the  stresses  and  strains. 
They  will  be  satisfied  by  a stress  function 2 F(x,y)  such  that 


(3) 


The  displacements  u(x,y)  and  v(x,y)  in  the  x-  and  y-directions 
respectively  are  necessarily  continuous  and  single-valued.  From  this 
condition  the  compatibility  equation  for  plane  stress. 


(U) 


^G.  B.  Airy,  Report  of  the  British  Association  for  the 
Advancement  of  Science,  1862. 
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is  derived. 

By  substituting  equations  (3)  into  equations  (l)  and  using 
equation  (U)  with  the  relation^ 
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the  differential  equation 


is  obtained  in  which^ 
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and 
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^March,  op,  cit.,  p.  7. 

^Far  wood  or  plywood  K,  as  defined  here,  is  probably  always 


greater  than  1 
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with 
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The  solution  of  (6)  is  facilitated  by  taking  F(x,y)  to  be  of 


the  form 


R[F,(y--+  i-«Ey)+Fx(xi-A.f£g?l 
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where 
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i is  the  imaginary  unit,  and  the  letter  R means  that  the  real  part  of 
the  expression  in  the  bracket  is  to  be  taken.  The  functions  F and 
are  any  analytic  functions  of  the  complex  variables  % -t  J.c\£  u and 


%+  J-  respectively . 
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CHAPTER  II 


GENERAL  EQUATIONS  FOR  STRESSES  AND  DISPLACEMENTS 

Two  types  of  problems  are  to  be  considered.  They  are: 

Type  I:  Problems  in  which  only  the  stresses  are  assigned 

on  the  boundary,  and 

Type  II:  Problems  in  which  certain  stresses  are  assigned 


on  the  boundary  together  with  certain  displacements 
on  the  boundary. 


In  both  types  the  load  will  be  symmetric  with  respect  to  the  y-axis  but 
shall  not  be  periodic.  Also  the  shear  stress  shall  vanish  on  the  entire 
boundary  in  both  types. 

Assume  then  a stress  function  in  the  form 


Differentiation  with  respect  to  the  proper  variables  yields 
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where 
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Noting  that 
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and  taking  the  real  part  of  the  expression  in  the  parentheses  makes 
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The  stress  Y is  symmetric  about  the  y-axis  and  in  order  to  satisfy 
this  condition  it  is  necessary  that 

A = fi  = o 

' ' 'Mil  ^12 

The  stress  function  then  reduces  to 
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The  stress  function  is  periodic  now  in  x and  in  the  form  of  a Fourier 
series.  The  stresses  are  to  be  non-periodic  in  x so  transform  the 
stress  function  into  a Fourier  integral,  that  is, 

oC  oO 
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where  A(t)  and  B(t)  are  to  be  determined  from  the  boundary  conditions. 
Differentiation  yields 
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The  stipulation  that  the  shear  vanish  on  the  boundary  makes 
it  necessary  that 
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Equations  (l6)  reduce  then  to 
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which  are  the  general  equations  for  the  stresses  in  the  form  of  a 
Fourier  integral. 

From  equations  (l) 
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Integrating  with  respect  to  x or  y as  the  case  may  be  yields 
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where  the  functions  u0(y)  and  vD(x)  are  determined  from 

bet  ()  <r-  _J_  y 

dij  Jv 

in  equations  (l).  These  are  the  general  equations  for  the  displacements 
u(x,y)  and  v(x,y)  in  the  form  of  a Fourier  integral. 
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CHAPTER  III 


SOLUTIONS  FOR  PROBLEMS  OF  TYPE  I 


Consider  now  problems  of  Type  I,  that  is,  problems  in  which 
only  the  stresses  are  assigned  on  the  boundary.  In  all  of  these 
problems  the  shear  stress  will  vanish  on  the  boundary  making 

X. r O UJ  he.  r\  Lj  - O 

which  will  permit  the  use  of  the  equations  obtained  in  the  preceding 
chapter. 

For  the  first  problem  assign  a normal  load  on  the  boundary 
distributed  in  the  manner  of 

H 

bUit* 

which  is  another  way  of  saying  that 
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since  the  distribution  of  the  stress  Yy  along  the  boundary  is  equal  to 
the  normal  load.  Then,  using  equations  (18),  when  y is  zero,  we  have 
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From  this  equation  it  follows  that 


or 
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Substituting  this  value  for  A(t)  in  equations  (18)  makes 
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Noting  that  the  integral  is  in  the  form  of 
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and  that  o( ^ = / , we  obtain  on  integrating  with  respect  to  t 
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In  order  to  integrate  these  last  integrals  notice  that 
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The  use  of  these  yield 
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which  are  the  stresses  produced  by  a load  distributed  in  the  manner  of 

Now  suppose  that  we  have  a point  load  of  magnitude  P applied 
at  the  origin  in  the  positive  y direction  instead  of  the  one  just 
discussed.  It  is  possible  to  get  the  stress  equations  for  the  point 
load  P from  the  stress  equations  just  obtained  by  allowing  the  constant 
b in  the  load  function  to  approach  zero  while  the  total  load  remains 
constant.  To  do  this  notice  that  the  total  load  is 
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and  let 
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Putting  this  value  in  the  stress  equations  (22)  and  letting  b approach 
zero  we  have 
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that  is,  the  total  force  on  a boundary  y=d  of  an  infinite  strip  in  the 
plate  mist  equal  the  force  P on  the  boundary  y=o  in  order  for  the 
strip  to  be  in  static  equilibrium.  The  shear  stress  vanishes  immediately 


when  yso  and 
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Let  us  also  investigate  to  see  if  these  stresses  reduce  to  the 
isotropic  case.  In  the  isotropic  case  o(.  ~€~/  so  noticing  that  when 

c^~ the  stresses  are  indeterminant  we  hold  constant  and  using 
1* Hospital's  rule  we  obtain 
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which  agrees  with  the  values  obtained  by  Timoshenko  if  the  choice  of 
axes  is  made  to  agreed  It  might  be  added  that  the  stresses  at  the 
origin  are  infinite  as  would  be  expected  from  a point  load  at  the  point 
of  application.  Therefore  equations  (2ii)  are  the  stress  equations  in 
a semi-infinite  plate  of  orthotropic  material  for  the  point  load  P. 

The  displacement  v due  to  the  point  load  P will  be  calculated 
now  since  it  is  to  be  used  in  the  next  chapter.  The  use  of  equations 
(l)  and  (2U)  yields 
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?S.  Timoshenko,  Theory  of  Elasticity,  New  York,  193U,  p.  8U. 
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The  function  v0(x)  must  either  be  a linear  function  of  x or  a constant. 
It  cannot  be  a linear  because  v is  symmetric  with  respect  to  the  y-axis . 
Therefore  v0  is  a constant  which  depends  on  an  assigned  deflection  at 
the  point  (0,d).  Then  the  equation  for  v due  to  a point  load  is 


where  k is  the  value  of  v at  the  point  (0,d). 


j ***** 

CHAPTER  IV 

fcf||||gt|g| 

SOLUTIONS  FOR  PROBLEMS  OF  TIPE  II 

Consider  now  problems  of  Type  II,  that  is,  problems  in  which 
certain  stresses  are  assigned  on  the  boundary  together  with  certain 
displacements  on  the  boundary.  Let  the  assigned  stresses  in  each  case  be 

V = 0 when  u-0 

Aj/  d (26) 

which  will  permit  the  use  of  the  equations  developed  in  the  second 
chapter.  Also  let  the  displacements  on  the  boundary  in  each  case  be 
the  displacement  v.  The  problem  shall  be  to  determine  the  load 
necessary  to  produce  each  of  these  boundary  conditions. 

Pbr  the  first  problem  assume  a displacement  along  the  boundary 

ir=  -tL 
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Then  from  equations  (19) 


and  when  y=  o we  have 
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The  function  v0(x)  must  be  zero  because  it  can  either  be  a linear 
function  of  x,  which  is  inpossible  due  to  the  symmetry  of  v,  or  a 
constant,  which  must  be  zero  so  that  v will  vanish  as  x becomes 
infinite.  Therefore 
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Using  equation  (2li)  we  have 
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r 


'ZTrr 


JL 


sjr^Ki 


I(&-  ^ 


(bt<*ey. 1 

(b+oizy)\Y>' 


(29) 


23 


When  y s o 
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which  satisfies  the  boundary  condition  equation  (27). 

In  order  to  find  the  load  distribution  consider  the  stress  ly 
where  from  equations  (18)  and  (28) 
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Using  equation  (22)  and  integrating  with  respect  to  t yields 
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In  order  to  integrate  this  last  integral  consider 


<-o 

& £-£<2.  (bvn. 


~k 

iiu+JL'-  ■ 


2k 


Differentiating  each  side  with  respect  to  k makes 
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We  can  use  this  result  to  obtain 
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This  is  the  distribution  of  the  stress  ly  along  the  boundary  and  there- 
fore the  distribution  of  the  load  which  is  a solution  to  the  problem. 

Consider  now  a semi-infinite  plate  whose  boundary  is  free  of 
shear  stress  and  that  under  the  load  the  displacement  in  the  y direction 
is  constant,  that  is. 


Xlj  ~ 0 U/hzn.  ^ = O 
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Then  the  effect  of  the  load  is  known  and  the  problan  is  to  determine 
what  load  will  produce  these  conditions  on  the  boundary  of  the  plate. 

In  order  to  determine  the  load  distribution,  first  suppose 
that  there  was  more  than  one  point  load  acting  at  various  points  on  the 
boundary  of  the  plate.  If  the  effect  of  each  load  taken  independently 
were  known  then  the  combined  effect  of  all  the  loads  would  be  the  vector 
sum  of  the  effect  due  to  each  load.  Keeping  this  in  mind,  consider  the 
total  load  as  a series  of  point  loads  which  are  distributed  from  a to 
-a.  The  point  load  will  vary  in  magnitude  so  that  the  combined  effect 
of  all  the  point  loads  will  produce  the  given  boundary  conditions.  Each 
point  load  will  be  a function  of  ^ , its  distance  from  the  origin 
along  the  x-axis.  Then  the  equation  for  the  displacement  v due  to  a 
point  load  a distance  f from  the  origin  is 
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by  using  equation  (25)  where  P(pd^  is  the  magnitude  of  the  point  load 
acting  at  the  point  ^,0.  The  constant  k in  equation  (25)  will  have  no 
effect  on  the  determination  of  P(^)  since  it  gives  a constant  deflection 
in  the  y-direction  for  every  point  along  the  boundary.  Hence  it  needs 
not  be  considered  here.  From  this  it  follows  that  the  equation  for  the 
displacement  v for  a series  of  loads  along  the  x-axis  from  a to  -a  would 
be  - a- 
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When  y =o 
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the  value  of  which  is  to  be  a given  constant  independent  of  x between  a 
and  -a. ^ Therefore  we  have  an  integral  equation  in  which  P(p  is  to  be 
determined. 

The  problem  is  similar  to  one  encountered  in  electrostatic 
potential  in  which  an  infinite  conducting  strip,  whose  width  is  2a,  is 
perpendicular  to  the  xy -plane. ^ If  the  strip  is  electrostatically 

^ It  is  easily  shown  that  a uniform  loading  of  the  straight 
boundary  from  a to  -a  would  not  produce  a constant  deflection  under  the 
load,  since  in  this  case  P(£)  would  be  a constant  and  equation  (33) 
readily  integrates  to  show  that  v varies  with  x when  y is  zero. 

70.  D.  Kellogg,  Foundations  of  Potential  Theory,  Berlin.  1929, 

p.  375. 
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charge^  then  the  potential  on  the  strip  mast  be  constant,  otherwise  a 
flow  of  charges  would  be  started.  It  is  shown  in  potential  theory  that 
the  potential  ((5  in  the  xy-plane  due  to  a long  straight  uniformily 
charged  wire  perpendicular  to  the  plane  is  logarithmic,  that  is, 

= 


< 


where  q is  the  charge  per  two  units  of  length  and  r is  the  radial 
distance  in  the  xy -plane  from  the  point  where  the  wire  intersects  the 
plane.  Then  think  of  the  strip  as  a series  of  wires  perpendicular  to 
the  xy-plane  along  the  x-axis  from  a to  -a.  Let  q be  the  total  charge 
per  two  units  of  length  of  the  strip  perpendicular  to  the  plane.  It 
is  then  necessary  to  determine  the  distribution  of  q along  x-axis  from 
a to  -a  so  that  the  potential  (f>  will  be  constant  along  that  part  of 
the  x-axis. 


If  we  use  the  complex  variable 

formation 


z - 


then  the  trans- 


maps the  charged  strip  on  the  x-axis  in  the  z-plane  into  the  unit 
circle  in  the  f -plane.  The  lower  half  of  the  strip  maps  into  the 
lower  half  of  the  circle  in  tte  f -plane,  the  upper  half  of  the  strip 
into  the  upper  half  of  the  circle.  The  rest  of  the  z-plane  can  be  mapped 
into  the  region  outside  of  the  unit  circle  in  the  f -plane.  Now  at 
great  distances  from  the  origin  the  potential  must  become  negatively 
infinite  like  -q  log  r.  At  great  distances  from  the  origin  f 
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becomes  infinite  like  - I ? / , that  is  like  — h • • Hence  the  function 

a.  1 1 a. 

<t>=-f  rlfl  wo 

which  is  constant  on  the  circle  | 91*1  # will  satisfy  the 

requirements  for  the  potential. 

It  is  also  shown  for  two-dimensional  potential  in  potential 
theory  that  on  the  charged  strip 


- H = -ln<r 

dK- 


(35) 


where  n is  the  normal  and  <T'  is  the  distribution  of  the  charge  q on 
the  surface  in  the  x-direction.  On  the  changed  strip  the  normal  is  in 
the  direction  of  y.  Therefore  equation  (35)  becomes 


Ai  - £j£ 
’ *1- 


(36) 


where  the  plus  and  minus  signs  mean  taking  the  partial  derivative  with 
respect  to  y from  the  positive  side  and  negative  side  respectively. 
Solving  for  j*  one  obtains 

Aji-lzl  (37) 


In  order  to  determine  whether  to  use  the  plus  or  minus  Sign  before  the 
radical  notice  that 


so  that  the  values  of  f will  not  be  inside  of  the  unit  circle  in  the 
f -plane,  for  no  points  in  the  z -plane  are  to  be  transformed  into 
points  inside  of  the  unit  circle  in  f -plane.  Then  when  investigating 
points  along  the  y axis  we  have  for  y =•  £ > o 


A 6 + 1 

j- 1 

" 

1 1 

A _£ 

<x 


+ [ 


S| 


which  shows  that,  for  this  last  quantity  to  be  equal  to  or  greater  than 
one,  the  plus  sign  must  be  chosen.  For  y=  £ < o we  have 


from  which  follows  the  negative  sign  is  to  be  used  when  y is  negative. 


Differentiation  yields 


which  on  the  strip  when  y is  zero  reduces  to 


Similarly  when  on  the  strip 
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r i - * . 

1 

1 

3T\> 

. V £ 2 /cC  - 

Then  using  equation  (35)  we  have 


2 / a.  - 


-ZnCT' 


or 

<r- 

^ ~ fa.'— 3cv  • 

This  is  the  desired  distribution  of  the  charge  q along  the  x-axis 
between  a and  -a. 


% 


Now  returning  to  the  original  problem  and  rewriting  equation 
(33)  in  the  form 

S3 


ir= 


illuZ  xvO 


2 n (f-«) 


££±\ P(  f )Jkj[(«ef*  f 


- <J_ 


2-  (/*-*)  J 


(38) 


P(f)  Aj  + 


a similarity  can  be  drawn  between  the  two  integrals  and  the  potential 
problem  just  discussed.  It  is  necessary  to  determine  the  load 
distribution  P(f)  so  that  the  integrals  will  be  constant  between  a and 
-a.  Also  the  functions  involved  are  logarithmic.  Then  assuming  that 
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P(p  is  distributed  in  the  same  manner  as  <T"  we  have 


P(f)  = 


K 

^ ”i  ^aF- V-1 


where 

-a. 

is  the  total  load  on  the  boundary. 

It  is  now  necessary  to  check  whether  or  not  this  does  satisfy 
the  boundary  condition  equation  (31).  Far  y = o equation  (38)  becomes 


a , 


or  remembering  that  o<  fl  = / 


~ir=  f A?  d-f 


‘ill"* 


- CL 


In  order  to  integrate  this  equation  let  f - Gl  £©o<p  > then  df  - 
and 
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The  second  integral  is  integrated  immediately,  but  in  order  to  find 
the  value  of  the  first  integral  consider  the  integral 

“ n 

where 


= I 


» 


Then  _ — 

__  — " I* 

This  integral  suggests  the  integration  of  the  complex  function 


.et* 


around  the  unit  circle  in  the  z-plane  since  on  the  unit  circle 

Z-  e/'e)  d-i  ~ d.» 

and 

= 2 to*3  & 

Rewriting  the  integral  we  have 


r 


J 

c 
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which  has  singularities  at 


= 0 And  2 = 


Now 


i- 


which  means  that  the  last  two  singularities  are  on  the  unit  circle. 
Each  of  these  will  be  excluded  from  the  contour  with  a circular  arc  of 
radius  £ and  with  the  center  at  the  singularity.  Now  noticing  that 
the  angle  O passed  through  during  the  integration  around  these 
singularities  is  the  same  for  each  one  since  the  radius  <5  is  the  sane 
for  each  and  letting 


we  have 


z-^r; 


Lirw  (z-Y'')  — ■■ 

1 i(i-r,xe-r, :)  r,(r,-n 


and 


3U 


Then  the  integral  of  f(z)  around  r^  and  ^ is 


-29 


I +-r; 


r,  ( rr  rJJ 


and  Z6t 


rJ  rr  n 


respectively.  The  sum  of  these  is 


which  reduces  to 


(rr-rJfcrt-O 


since 


Therefore  using 


it' 


Li  171  ^ -2L(l  + 

2--*0  Z(-kz"-2iz-i-'k) 


-JL 

'k 


we  have 
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or 


n 


— 2 $ 
I - 


The  is  an  even  function  so  this  reduces  to 


Integrating  this  with  respect  to  k in  the  following  manner  we  have 


zA£^t- sL<f  = I *1 

I 6 I 


Wien  k =^1 


which  reduces  to 

77 


J 

0 
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Let  4-lX.  d.<f-  2JLI 


then  we  have 


„ r % 

•?  77  9-  + J /m™.  /Ldt 


or 


2 Ti  1 - H T\  2 - ~2n  ^ 


Hence 

({~iz  CdnCffcitf  - 2 ,7  Mg  2 -k 
o 

and  equation  (39)  becomes 


-ir  = 

^ T,  **€ 


or 


"IT 


= u.  ju. 


^7T'£ 


Therefore  the  boundary  condition  equation  (31)  is  satisfied, 
the  load  is  the  negative  of  the  stress  on  the  boundary  and,  since  there 
is  no  load  for  , equation  (32)  is  satisfied,  and  the  original 

stress  function  was  determined  so  that  there  was  no  shear  stress  on  the 
boundary,  satisfying  equation  (30).  The  assumed  load  is  therefore  a 
solution  to  the  problem. 
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